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Abstract. We develop a local theory of lacunary Dirichlet series of the form 

oo 

C;. exp(— 2(/(fc)), Re(2) > as z approaches the boundary iM, under the 

k = l 

assumption — > oo and further assumptions on Cj. . These series occur in many 
applications in Fourier analysis, infinite order differential operators, number 
theory and holomorphic dynamics among others. For relatively general series 
with cj, = 1, the case we primarily focus on, we obtain blow up rates in measure 
along the imaginary line and asymptotic information at z = 0. 

When sufficient analyticity information on g exists, we obtain Borel sum- 
mable expansions at points on the boundary, giving exact local description. 
Borel summability of the expansions provides property-preserving extensions 
beyond the barrier. 

The singular behavior has remarkable universality and self-similarity fea- 
tures. If g{k) = k^, Ci. = 1, b = n or b = (n + l)/n, n G N, behavior near the 
boundary is roughly of the standard form Re(2)~' Q(x) where Q(x) = 1/g if 
X = p/q € Q and zero otherwise. 

The Botcher map at infinity of polynomial iterations of the form Xn+i = 
XP{xn), |A| < Xo{P), turns out to have uniformly convergent Fourier expan- 
sions in terms of simple lacunary series. For the quadratic map P{x) = x — x^, 
Ao = 1, and the Julia set is the graph of this Fourier expansion in the main 
cardioid of the Mandelbrot set. 
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1. Introduction 



Natural boundaries (NBs) occur frequently in many applications of analysis, in 
the theory of Fourier series, in holomorphic dynamics (see [8], [7], [5] and references 
there), in analytic number theory, see [20], physics, see [13] and even in relatively 
simple ODEs such as the Chazy equation [llj . an equation arising in conformal 
mapping theory, or the Jacobi equation. 

The intimate structure of NBs turns out to be particularly rich, bridging analysis, 
number theory and complex dynamics. 

Nonetheless (cf. also [S]), the study of NBs of concrete functions is yet to be 
completed from a pure analytic point of view. The aim of the present paper is a 
detailed study near the analyticity boundary of a prototypical functions exhibiting 
this singularity structure, classes of lacunary series. For such functions, we develop 
a theory of generalized local asymptotic expansions at NBs, and explore their con- 
sequences and applications. The expansions are asymptotic in the sense that they 
become increasingly accurate as the singular curve is approach, and in many cases 
exact, in that the function can be recovered from these expansions. 

Lacunary series, sums of the form h{s) = '^j>iCjS^^ , or written as Dirichlet 
series f{z) = J2j>i Cje~^^^^^), where j / gj = o(l) for large j, often occur in applica- 
tions and have deep connections with infinite order differential operators, as found 
and studied by Kawai fT2], [13j . Under the lacunarity assumption above, if the unit 
disk is the maximal disk of analyticity of h, then the unit circle is its NB ([15j). 
For instance, the series 



studied by Jacobi [10] before the advent of modern Complex Analysis, clearly has 
the unit disk as a singular curve: h{s) — > +oo as \s\ t 1 along any ray of angle 
2^"r7i7r with m,n eN. 

We show that if Cj = 1 and g' — > oo, then, in a measure theoretic sense, f(x + iy) 
blows up as X ^ 0^ uniformly in y at a calculable rate. We find interesting 
universality properties in the blow-up profile. 

In special cases of interest, Borel summable power series, in powers of the dis- 
tance to the boundary, and more generally convergent expansions as series of small 
exponentials multiplying Borel summed series power serie|3 representations can be 
determined on a dense set on the singularity barrier. Examples are 



These are simple instances of Borel (or Ecalle-Borel) summed transseries. A brief summary 
of the definitions and properties of transseries and generalized Borel summability, and references 
to the literature are given in i|5.6l 



oo 




(1) 
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• Sj>i where 6 > 1, or its dual g > 1, where b^^ + q^^ — 1, is integer 

(relating to exponential sums and van der Corput dualities |17j : the special 
self-dual case J2j>i ' related to the Jacobi theta function); 

More generally, if Cj — c{j) and gj — g{j) have suitable analyticity properties in 
j, then the behavior at z — 0^ , and possibly at other points, is described in terms 
of Ecalle-Borel summed expansions (see footnote [T]) . Then the analysis leads to a 
natural, properties-preserving, continuation formulas across the boundary. 

In general, the blow-up profile along the barrier is closely related to exponential 
sums, expressions of the form 

N 

Sn^YI ^» exp(27ri5(n)), g{n) G M (2) 

where for us 5' — > cxs as n ^ oo. The corresponding lacunary series are in a sense 
the continuation of ^ in the complex domain, replacing 27ri by —z, Re(2:) > 0, and 
letting N ^ 00. The asymptotic behavior of lacunary series as the imaginary line 
is approached in nearly-tangential directions is described by dual, van der Corput- 
like, expansions. The method we use extend to exponential sums, which will be the 
subject of a different paper. 

2. Results 
2.1. Results under general assumptions. 

2.1.1. Blow-up on a full measure set. We consider lacunary Dirichlet series of the 
form 

00 

/(z) = ^e-^« (3) 

but, as it can be seen from the proofs, the analysis extends easily to series of the 
form 

00 

/(2) = ^c,e-^^« 

under suitable smoothness and growth conditions, see ^3.8\ The results in the paper 
apply under the further restriction. 

Assumption 1. The function g is difjerentiable and g'{j) 00 as j ~f 00. 

In particular, g is eventually increasing. By subtracting a finite sum of terms 
from / (a finite sum is clearly entire), we arrange that g is increasing. If g(0) = a, 
we can multiply / by e~^° to arrange that g{0) — 0. 

Normalization 1. (i) g is differentiable on [0,oo), g' > and g' ^ 00 along M+. 
(n) 5(0) = 0. 

Notation. We write \H{-) \ = l+o(l) if \H{y)\dy converges to the Lebesgue measure 
dn{y). 

Under Assumption [U after normalization, we have the following result, giving 
exact blow-up rates in measure, as well as sharp pointwise blow-up upper bounds. 
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Theorem 1. (i) We have the uniform blow-up rate in measur^. 

/"OO 

+ / e-^s^'>ds{l + o{l)) (4) 



It can be checked that e-^3(s)x^g > 3-1(1/2;) 

— > (X) as X ^ 0+; see also Note[l\ 

(ii) The following pointwise estimate holds: 

poo 

\\f(x + f)\\oo< e-ff(^)^ds(l + o(l)) asx^O+ (5) 



This is sharp at z — 0, cf. Proposition and in many cases it is only reached at 
z = 0; see Proposition^^ 

2.1.2. General behavior near z ^ 0. At z = 0"*" a more detailed asymptotic descrip- 
tion is possible. 

Theorem 2. (i) As z ^ O'^ we have 

/(z)= / e-^^^'^'ds - - + 0(1) ^00 as 2^0+ (6) 

Jo 2 

In fact, 

f{z)- e-^^^'Us^-z e-^"{,g-i(u)}dw (7) 



(where {•} denotes the fractional part, and we used g{0) — oj^. 

(ii) If g{s) has a differ entiable asymptotic expansion as s 00 in terms of 
(integer or noninteger) powers of s and logs, and g{s) ^ const. s^ , b > 1, then after 
subtracting the blowing up term, f has a Taylor series at z = (generally divergent, 
even when g is analytic, which can be calculated explicitly), 

f{z)= e-'s^''>ds- -+zs{z), seC°°[0,oo) 

(as an example, see h61\) . 

Note 1. Often g has an asymptotic expansion starting with a combination of pow- 
ers, exponentials and logs. Let (j) = g~^ . Then 4>(vx) / 4>(v) 0i(a;) as v ^ 00 and 
X > Q is fixed and 

00 poo 

e"'^^'Us^Cgg-^{l/x){l + o{l)); Cg^ e-''(t>i{u)du 

Jo 

For instance, if g{k) = fc'', 6 > 1 we have, as p — > 1, 

\f{z)\=x-^^ r(i + i/fe)-i/22-A(i + o(i)) (.T^0+) 

/(x)=x-i r(l + l/6)(l + o(l)) (x^0+) (8) 



^It turns out that in general, or a.e. convergence of |/| do not hold. 

■^As mentioned, often this maximal growth is achieved at zero but in special cases it occurs, 
up to a bounded function, on a dense set of measure zero. 
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2.1.3. Blow-up profile along barrier. 

Theorem 3. (i) Assume that for some y G R there is a smooth increasing function 
p{N]y) =: p{N) G (0, iV] such that the following weighted exponential sum (see 
|17j ) has a limit: 

N 

Sp^N := p{Ny'^ ^'^^^^'^ ^(2^) as N^oo (9) 

where p" is uniformly bounded and nonpositive for sufficiently large k. (Without 
loss of generality, we may assume p"{k) < for all k.) Let 

$(x) ^ / e"^ff("V'Hrfw (10) 



Then, we have the asymptotic behavior 

fix + iy)^L{y)^x)+o{Mx)) (x 0+) (11) 
(a) As a pointwise upper bound we have: 

liiasnp \S p^n\= L <oo ^ fix + iy) ^Omx)) {x ^ 0+) (12) 

JV>0 

2.2. Results in specific cases. 

2.2 A. Settings leading to convergent expansions. The cases g{j) = and g{j^ = a^, 
a > 1 are distinguished, since the expansions at some points near the boundary 
converge. 

Proposition 1. Ifb— 2, then we have the identity 



* * fc=i 

Clearly, this is most useful when z ^ 0. It also shows the identity associated to 
the Jacobi theta function 

^^■'^-2i--2^rj[-^) (1^) 

Proposition 2. If g{j) = , a > 1, then, as z 0^, f{z) is convergently given 
by 

^ "bi^ + ^ n!(l - a") + + bi^ ^ 

= J^ + f2-t^ + co-^ r log^[-is/C)^]e-^ds (14) 
loga ^-^ n\(l — a"-) 2TTt In 

where C. = 1 — e^^. (Here logjj is the usual branch of the log with a cut along 
and not the log on the universal covering ofC\ {0}.) 

It is clear that for a G N the transseries can be easily calculated for any z = 
pexp{2'!Tim/a^), {m, j) G N, < p < 1 since 

3 

f{peyip{2mm/a^)) = ^ exp(27r«ma"^^) + f{p) 
where the sum is a polynomial, thus analytic. 
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2.2.2. Borel summable transseries representations; resurgence (cf. §5. 6\} . When 
Cj = 1 it is clear that the growth rate as z — > 0+ majorizes the rate at any point on 
zM. There may be no other point with this growth, as is the case when g{j) — , 
b E (1,2) as seen in see Proposition [5] below, or densely many if, for instance, 
g{j) — , 6 G N, or when g{j) — , a E N. The behavior near points of maximal 
growth merits special attention. 

For g = j^, 1 < b 2, f has asymptotic expansions which do not, in general 
converge. They are however generalized Borel summable. 

Define d by 

1 + 5^' 

Theorem 4. Let g — j^; b > 1. Then, 

(i) The asymptotic series of f{z) for small z, 



1 \ 1 . 

n 1 . 



/o = r ( 1 + - J z-t - - + — - {-iy'')ajb + l)b,z^ (16) 



is Borel summable in X ^ z ^Z*-^ ^\ along any ray a,Tg{X) = c if c ^ ~ argk'^s ±, k E 
N, where s± =t±^ 2Tiit'"£' and t± = (±27rj/6)''/(^"^) . More precisely, (a) 

1 

/(z) -r(l + 1/6)^-1/" - i + / e-'-''''''"''H{s)ds (17) 

2 Jo 

where H{s) — Hi,{s^^^), where Hi, is analytic at zero and -fffc(O) = 0; (b) H is 
analytic on the Riemann surface of the log, with square root branch points at all 
points of the form k''s±,k E N and (c) making appropriate cuts (or working on 
Riemann surfaces), u~^^~^^ /^H is bounded at infinity. 

//arg(X) E (0_,0+), then CBfo = / /n a general complex direction, f has a 
nontrivial transseries, see (Hi). 

(Hi) For a given direction ip, a be ±1 if ±ip > and otherwise. If aargz E 
{9a-, 7^/2), then the transseries of f is 

oo oo 

f^(z) + a—Y,e~"' ^Cj-(afc) ^("-d z^^ (18) 

k=l j=o 

and it is Borel summable as well. 

Note 2. The duality k^ ^ k'^ is the same as in van der Corput formulas; see [17]. 
2.2.3. Examples, (i) For b — 3, f the transseries is given by 



'!Tik \ " 1 5 (ik) 4 1 

Z~3 H 3 5 ^^ + 

6 J 32 6 4 TT 4 



(19) 



where with = 7r'^/^y^||(— 1) 4 and 

r(4/3) 1 z 



fo{z) 



ji/3 2 120 792 



Note that the variable of Borel summation, or critical time, is not 1/z but 2" 



-1/(6-1) 
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(ii) For b — 3/2, with k — 32i7r'^/27, / has the Borel summable transseries 

oo 

/o(z) + (7^e' 
fc=i 

where 9± = 7r/4 and 

2.2.4. Properties-preserving extensions beyond the barrier. It is natural to require 
for an extension beyond the barrier that it has the foUowing properties: 

(i) It reduces to usual analytic continuation when the latter exists. 

(ii) It commutes with all properties with which analytic continuation is com- 
patible (principle of preservation of relations, a vaguely stated concept; this 
requirement is rather open-ended). 

Borel summable series (more generally transseries) or suitable convergent represen- 
tation representations allow for extension beyond the barrier, as follows. In the 
case g{j) = (and, in fact in others in which g has a convergent or summable 
expansion at infinity), f{z) can be written, after Borel summation in the form (see 

'Hi{s)ds ^ zT^-'^ e-*i/i(tzi/(^-i))dt=: / e-''F{s,z)ds 
Jo Jo 

(21) 

where Hi is analytic near the origin and in C except for arrays of isolated singulari- 
ties along finitely many rays. Furthermore, Hi is polynomially bounded at infinity. 
This means that the formal series is summable in all but finitely many directions 
in z. 

Definition 3. We define the Borel sum continuation of f through point z on the 
natural boundary, in the direction d, to be the Borel sum of the formal series of f 
at z in the direction d, if the Borel sum exists. 

This extension simply amounts to analytically continuing F{s, z) in z, in h21]) . 
for small s and then analytically continuing F{s,z) in s for fixed z. along M^. 

2.2.5. Notes, (cf. Appendix gS^l) 

(i) The Borel sum provides a natural, properties-preserving, extension [4j. Borel 
summation commutes with all common operations such as addition, multi- 
plication, differentiation. Thus, the function and its extensions will have the 
same properties. 

(ii) Also, when a series converges, the Borel sum coincides with the usual sum. 
Thus, when analytic continuation exists, it coincides with the extension. 

(iii) With or without a boundary, the Borel sum of a divergent series changes 
as the direction of summation crosses the Stokes directions in C. Yet, the 
properties of the family of functions thus obtained are preserved. There may 
then exist extensions along the barrier as well. Of course, all this cannot 
mean that there is analytic continuation across/along the boundary. 

See also Eq. (|14p . a convergent expansion, where z can also be replaced by — z. 
In this case, due to strong lacunarity, the extension changes along every direction, 
as though there existed densely many Stokes lines. 



4\/2((Ti) 2 7rfc2 (o'i)2 z^ 
3 z 16V27rt fci 



(20) 
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Figure 1. The standard function Q{x). The points above Q 
0.05 are the only ones present in the actual graph. 




e.2 



Figure 2. Point-plot graph of (34,4, normalized to one. 



2.3. Universal behavior near boundary in specific Ccises. In many cases, 
<t){g~^{u/ x)) has an asymptotic expansion, and then $ in turn has an asymptotic 
expansion in x. Detailed behavior along the boundary can be obtained in special 
cases such as g{j) = j^; 6 e N, or g{j) = 6 € N. Properly scaled sums 

converge to everywhere discontinuous functions. 



Theorem 5. (i) If g{j) = j^, b eN, we have limsupa;^/''|/(a; + iy)\ < M < oo. 
We let Qb Js) = lim x^^'^fix + 2ms) (whenever it exists); then, 

x->0+ 
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•25 



Figure 3. Point- plot graph of Qs^s; Q3,2 follows from it through 
the transformation (1241). 



Sb,b 



1 " 

-E 

n ^- — ' 



r(i 



a.e. 



m 

e n , s = — , m,n £ N relatively prime 
n 



(22) 



(ii) For 6 = 3/2 and z — > / grows like z ^Z'"^. For any other point on iR the 
growth is slower, at most (Rez)^^/^. Furthermore, 



Q3/2,l/2{Vs) 



1 



E 

1=1 



m 

-2TTi — 

e n 



16 n 
27to 



, n,m as in (i) 



(23) 



a.e. 

/n particular, we have the profile duality relation 



r(4/3) """^^ V- " " y (24) 

for any s G M /or which (33/2(5) and/or Q3,3{s) is well defined, for instance, in the 
cases given in and h2S\l . 

For large n, the sum over I in Eq. is, statistically, expected to be of order 
one. After x"^ rescaling, the template behavior "in the bulk" is given by the familiar 
function Q{x) = if x = m/n, {m,n) G and zero otherwise, shown in Fig. [31 

2.4. Fourier series of the Botcher map and structure of Julia sets. We 

show how lacunary series are building blocks for fractal structures appearing in 
holomorphic dynamics (of the vast literature on holomorphic dynamics we refer 
here in particular to [5], [7] and [16j). In ^we mention a few known facts about 
polynomial maps. Consider for simplicity the quadratic map 



(25) 
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It will be apparent from the proof that the results and method extend easily to 
polynomial iterations of the form 

Xn+l = \Pk{Xn) (26) 

with A relatively small. The substitution x ~ — (Ay)~^ transforms (|25p into 

A(l + Vn) 

By Botcher's theorem (we give a self contained proof in 21 for ([25]) . which 
extends in fact to the general case), there exists a map analytic near zero, with 
0(0) = 0, 0'(O) = A"^ so that ((^/(/)"^)(z) = z^. Its inverse, i/', conjugates ^} to 
the canonical map z„+i = z^, and it can be checked that 

^{zf = A^(z2)(i + v-(z)); ^(0) = 0, ^'{Q) = A (28) 

Let ^(B) denote the Banach space of functions analytic in the unit disk D and 
continuous in D, with the sup norm. We define the linear operator X = T2, on 
Am by 

-J 00 

(T/)(z) = -^2-V(^''') (29) 

This is the inverse of the operator / ^ 2/ — J^^, where f^'^{z) = f{z^)- Clearly, Tf 
is an isometry on ^(D) and it maps simple functions, such as generic polynomials, 
to functions having 9D as a natural boundary; it reproduces / across vanishingly 
small scales. 

Theorem 6. (i) tfj and H = l/i/j are analytic m fA, m D x D ('A G D corresponds 
to the main cardioid in the Mandelbrot set, see and continuous m D x D. The 
series 

00 

V'(A,z) = Az + ^zA'=+Vfe(2) (30) 

converges m D x D (and so does the series of H), but not m D x D. Here 

00 

V;i(z) = Tz = ^2-'=+iz2'= (31) 

k=l 

(note that ip{z) = 5^ + 5 Jq s^^h{s)ds with h given in U^), and in general 

(k-l k-1 \ 

z ^ ipj'^ipk-i-j - ^ V'j Vfe-j (32) 
3=0 i=i / 

(ii) All xj}k,k > 1 have binarily lacunary series: In ipk, the coefficient of zP is 
nonzero only if p has at most k binary- digits equal to 1, i.e., 

p = 2^1 H h 2^" , ji = 0, 1, 2, ... ( "k "is the same as in ipk) (33) 

(Hi) For |A| < 1, the Julia curve of h25\} is given by a uniformly convergent Fourier 
series, by (i), 

J = {- (Rei7(e"),ImiJ(e'*)) :i G [0,27r)} (34) 



^The notations tp, H, designate different objects than tiiose in previous section. 
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Figure 4. The Julia set of Xn+i = Aa;„(l - x„), for A = 0.3 
and A = 0.3«, calculated from the Fourier series ((30|) discarding all 
o(A^) terms. They coincide, within plot precision, with numerically 
calculated ones using standard iteration of maps algorithms. 

Remark 3. The effective lacunarity of the Fourier series makes calculations of the 
Julia set numerically effective if A is not too large. 

Note 4. 

(i) Lacunarity of ipk is a strong indication that ip has a natural boundary (in 
fact, it does have one), but not a proof. Transseriation at sm(7w/ar points and 
summation of convergent series (here in the parameter A) do not commute. 
The transseries of H on the barrier can be calculated by transasymptotic 
matching, see but in this case is more simply found directly from the 
functional relation; see Note [71 

(ii) In assessing the fine structure of the fractal using the Fourier expansion trun- 
cated to o(A"), the scale of analysis cannot evidently go below 0(A"). 

Remark 5. (i) For |A| sufficiently small, Theorem [6] provides a convenient way 
to determine the Julia set as well as the discrete evolution on the boundary. 

(ii) For small A, the self similar structure is seen in 

^,(pexp(2W2")) = E P"^^^^-) + gi^p) 

k=l 

where the sum is a polynomial, thus analytic. Up to a scale factor of 2""+-'^, 
if |A| <C 2~"+^, the nontrivial structure of ip at exp(27riTO/2") and at 1 are 
the same, see NotelU that is 

V'(exp(27rim/2")) = 2""+ V(l) + regular + o(A) 

Exact transseries can be obtained for ip; see also Note[71 

(iii) For iterations of the form Xn+i = X'^~^Pk{x) where Pk is a polynomial of 
degree k > 2, the calculations and the results, for small A, are essentially the 
same. The lacunary series would involve the powers z^^ . For instance if the 
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recurrence is Xn+i = A^a;„(l — 3xn + x'^), then ^ is to be replaced by the 
solution of 

and the small A series will now have ipi — X^fc^i 3~'^~"'^z'^'' and so on. 



3. Proofs 



Lemma 4. Under the assumptions of TheoremUl (Hi), h{y) := g ^(y) also has a 
dijferentiable asymptotic power series as y oo. 



Proof. Straightforward inversion of power series asymptotics, cf SJS] 



□ 



3.1. Proof of Theorem [TJ The proof of (|4]) (i) essentially amounts to showing 
that l/p is diagonally dominant, in that terms containing g{j) and g{k) with j ^ k 
are comparatively small, as shown in ^ 33. 81 

(ii) Since g is increasing on (0,oo), the result follows from the usual integral 
upper and lower bounds for a sum. Equation ^ follows from simple calculations, 
cf. M 



3.2. Proof of Theorem[2l We prove part (ii); part (i) is similar, and simpler. By 
standard Fourier analysis we get 



1 sin 2k'KU 



where 



1 



k=l 



sin 2kTTU 



k=l 



kiT 



kn 



< C 



M 



Z 



Si(7r) as M ^ oo 



(35) 



(36) 



where Si is the sine integral, and Si(7r) is the constant in the Gibbs phenomenon. 
Let gm be an analytic function such that g(s) — gm{s) = o(s^™) for large s. 

Lemma 5. The analysis reduces to the case where h is a finite sum of powers. 
Indeed, 



f 



+ / e-''"g-\u)du + 



'{hM}du + Rrn-l{x) (37) 



9(0) 



where Rm-i is C™ ^ and is a truncation of the asymptotic expansion of h, 
such that h{u) — hrn(u) — o{u~™). 

Proof. We have 



/ e-^"{/i(u)}dM = V / e-^"(/i(M) - N)du 

= E/ e--^{hUu)~N)du + R,^_,{x) (38) 
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where 

rgmiN) /■3(A'+l) 



R,n-i{x) = V / + / e-^^{h{u)}du 

+ Y. e-''''{h{u)-hra{u))du (39) 

Using jSSl) and ^ the proof foUows and the fact that g{N)-gm{N) = o(iV~'") and 
h{u) — h.m{u) = o(u^™), the sum is rapidly convergent, and the result follows. □ 

Lemma 6. If h is a finite sum of powers, then f — e~^^(°) — e'^^^g^^{u)du £ 
C°°. 



Proof. Using ([35]) and ([36]) we have 

°° j,„sin2fc7r/im(u) 

A.— 1 



nOO 1 poo 

fi{x) = X / e-^«{/i™(^.)}du = :r " V / e- 

^-Ej^E-i '--v. (40) 

We deform the contours of integration along the directions aa respectively, say for 
a — tt/2. The integral 



/i ^ / u"e2'='^'''-(")du (41) 
Jo 

exists for any n and it is estimated by 

< const.k-''^''+^^T{n{b + 1)) (42) 







The termwise nth derivatives at 0^ of the series of fi{x) converge rapidly, and the 
result follows. □ 

3.3. Proof of Theorem [3l 

Proof. Let e > be arbitrary, let 

N 

i=i 

and let Ni be large enough so that \SN/piN) - L\ < e for N > Ni. Then, by 
looking at e**^/ if necessary, we can assume that L > 0. We have 

oc oo ^ 

f{x + *2/) = V e-^a^'^\Sk ~ Sk-i) = V -fr (e-^fW - e-^^^'^+^^)p{k) 
k=i P^''^ 

OO OC 

= L{y)J2ie'''^^'''^^ - e-^f('=))p(fc) + ^(e-^s^'^) - e-''s^''+^^)dkp{k) (43) 

fc=l k=l 

where dk — L{y) — Sk/p{k) as fc — > oo. Now, 



5]e--sW(p(fc + l)-p(fc))^ / 
k=i J° 



-''3('''^p'{k)dk $p(x). (44) 
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and under the given assumptions $p — )■ cx3 as x ^ 0. Note that $p ^ cx) as x — > 
and the facts that e'^'a^'''' - e-^fC^+i) > 0, p{k) > and dk 0, readily imply that 
the last sum in (|43|) is o($p(a;)); p2p is follows in a similar way. □ 

3.4. Proof of TheorerE(5j We rely on Theorem[3l and analyze the case b = 3/2; 

m 

the case 6 £ N is simpler. We have (3 — t/(27r). Let Sj = ^ e jt- . It is clear 

k=l 

m 

that for m,n € N we have j~ Sj ^ Yl, ^ ^ = ^mn- On the other hand, by 

fc=i 

summation by parts we get 

3p few - 1 3/i ^ 

and it is easy to check that, as — > oo we have 

I 3/3 (1 + 1/3) I few imn = YT ^ 

where we used the definition of fcjv following eq. (|45)) . which implies /cyv ^ |A^^/^/3. 
The result follows by changes of variables, using pTjl and noting that 

It is also clear that Ij^^^jl < 1 and a similar calculation provides an overall 
upper bound. 

Section 15.41 provides an independent way to calculate the behavior along the 
boundary. 

3.5. Proof of Proposition [TJ For b — 2 using (|40l) and e-P"" sin 2knp 2 dp = 
ic'^'^ke^^^^ x"^!'^ we immediately obtain 



TT 1 VTT'^ 

e ^ 



2~^x 2 V 2^ 

3.6. The case gij) = cu' • It is easy to see that, for this choice of we have the 
functional relation 

f{z)- f{az) = e-^' 

Since ^ 1 as z — ^ 0+, the leading behavior formally satisfies f{z)~f{az) ^ 1 

i.e. / z ~ 

log a 

log z 

In view of this we let f(z) — h G(z) which gives G(z) — Giaz) = — 1 

log a 



G(z) = G(z/a) + (l-e-"/'') 



(48) 
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We first obtain a solution a solution G of the homogeneous equation and then write 
G = G + h{a^), where h now satisfies 

h{y)^h{y + l) (49) 

Iterating (^5]) we obtain 



oo oc 



" ' - n— 1 ^ ^ 



n=0 k=l 



which is indeed an entire function and satisfies the functional relation. Now we 
return to (|49p . Since by its connection to f{z), h is obviously smooth in y, it can 
be expressed in terms of its Fourier series 



A;— — oo 

where the coefficients Ck — C-k can be found by using the original function /. Recall 
that we have 

f[z) = + G{z) + h 

log a \ log a 

which implies 

hiy) = f{ay)+y-G{ay) 
Lemma 7. The Fourier coefficients of the periodic function h are given by 

c^J-rf-p^) (MO) (50) 



log a \ log a 



co= / f{ay)dy+ f ydy~ f G{ay)dy 
Jo Jo Jo 



/ 1 -lkl„2 

Since Ck ^ \ — e '"K" , (implying that the Fourier expansion for h is valid 

V Trloga 

exactly for Re (z) > Oj we have 

^ = " 1^ + S ^^!(T^ + '^^ + gH^i^J 

The series further resums to 

-1 oo / \n poo 

/(^) = -p^ + E ,7 n, +-^-i^ / logK (-.'^) e-d. (51) 
logo rt! (1 — a") zm Jo ^ ^ 

valid for z > 0. 

The proof is given in 
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Similar results hold for other rational angles if 6 G N (by grouping terms with 
the same phase). For example, 

oc 

> g-2"(z+f77i) 



E 



oo oo oo oo 



e 

n— n— n— n— 



oc 



n— n— n— n.=0 

(52) 

3.7. Proof of Theorem |4l (i) and (ii): By an argument similar to the one leading 
to ([40]) we have 

/>oc -1 -1 paioc 

7o 2 ^^2fc^z^f^^7o 



The exponential term ensures convergence in k. Taking = t we see that 

e-'"'''dp = r{l + l/b)z-^^'' (54) 
We now analyze the case cr = — 1, the other case being similar. A term in the sum 



IS 



1 r ^-.„-2fe.»V.^^^ A\ /-^^.,(t+2.,tV.)^^ 



2kTri ./n V z 

d 



1 

zj Jo 



' pOD poo 

/ e-'"'''H^{s)ds^ z-'^ e-^"'^"'^""7I„(s/fc'')ds (55) 
Jo Jo 

where Vk = k'^z-^/^^^^\ H_{s) = (t>'{s) and 

^{(t){u)):^(t>{u) + 2m(j){uf'^ =u (56) 
(a) Near the origin, we write (/> — u^Hi,, and we get 

ZTTt 

and analyticity of Hi, in u'^^^ follows, for instance, from the contractive mapping 
principle. 

(b) The only singularities of 0, thus of are the points implicitly given by 
0. 

(b) It is also easy to check that for some C > we have 

uniformly in C, with a cut at the singularity, or on a corresponding Ricmann surface. 
From ([57]) we see that the sum 

oc 

if_ = yij_(./fc'') (58) 

fc=l 
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converges, on compact sets in s, at least as fast as const ^ fc^**, thus it is an analytic 
function wherever all H_ [s/k'^) are analytic, that is, in C except for the points k'^s^. 
Using an integral estimate, we get the global bound 

\H\ < consty , „ , < const|u|(^-i)'/^ (59) 

as |u| —-^ oo. 

The function H in the lemma is simply iJ_ + iJ_ where Hj^ is obtained from 
by replacing i by —i. The calculation of the explicit power series is straightforward, 
from (|56p , (|58p and the similar formulae for , using dominated convergence based 
on We provide the details for convenience. 

We write the last term in in the form 

b 

U{z)^(^^''' je-^'^^^H.ds (60) 
c 

where the contour C is a curve from the origin to oo in the first quadrant. 
Watson's lemma implies 

N 

fk{z) = Y,hjk-^^7J-^ + RN{k,z) 

where 6jCan be calculated explicitly from Lagrange-Biirmann inversion formula used 
for the inverse function (j)^^ and RN{k, z) ^ Ck~^^^^^^ z'^ , for arbitrary G N. It 
follows that 

j{z) ^ r(i + -)z"i _ - + A _ {-iy'^)ajb + i)h,z' (61) 

i=i 

which holds for z — s- in the right half plane in any direction not tangential to the 
imaginary axis. 

(iii) We obtain the transseries (which gives us information near the imaginary 
axis) by using the global properties of H- , and standard deformation of the Laplace 
contour. 

As z goes around the complex plane, as usual in Laplace- like integrals, we rotate s 
in (j83p simultaneously, to keep the exponent real and positive. In the process, as we 
cross singularities, we collect a contribution to the integral from the point s_ above; 
the contribution is an integral around a cut originating at s_. The singularity is 
integrable, and collapsing the contour to the cut itself, we get a contribution again 
in the form of a Laplace transform. This is the Borel sum (in the same variable, 

Generically s_ is a square root branch point and we have 

2^) 



ds- \ - _ 



The asymptotic expansion of the cut contribution, Borel summable as we men- 
tioned, is 



exp -(2^)^(6T^-6T^)zT^ f-j E^Hf)""" 
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The exponential term ensures convergence in k of the Borel summed transseries. A 
similar result can be obtained for —k. 
Thus, the transseries of / is of the form 

fiz) = r(i + - 1 + ^ f;(i - i-iy')ajb + mz^ 



+ < 



°° , N / 2» s-ir^ / fr' MT^T- °^ -236 + 2-6 23-1 

^ ^ e(''-i)(Tf)''-MV)'-^ J2 cjk ^c-^) z^^ -^^avgz^Oi 



±. £ e('>-^m)^i^)^ £ e,(-fc)-W.« e, < arg. < 

fe=l ]=0 

Remark 6. The analysis can be extended to the case 

00 

f{z) = J2me-'''' (A^O) 

k=l 

by noticing that 



POO 

fiz)= / izF{p)-F'{p))e-''^\pi]dp 
Jo 



If the expression of F{k) is simple, for example F{k) is a finite combination of terms 

of the form fj,''^ fc'^(log fc)"* (m = 0, 1, 2, 3...), the method in this section applies with 
little change to calculate the transseries of f{z). 

For special values of b, asymptotic information as z approaches the imaginary 
line can be obtained in the following way. Let z = 5 + 2iTif3 and 



f{z) = ^e-^'^^+^-^z') (Re(^) > 0) 



fe=i 

m 

If 1 < 6 G N, we may obtain the asymptotic behavior for all rational (3 = — , by 

n 

noting that e~^ (2m0) _ g-(fe+n) (27rj/3) ^^^^ splitting the sum into 

/(-) = E E e"^"'^'^ - ) = E E ^-^"^'^'^^^ 

j=o 1=1 1=1 3=0 

It follows that 



^ fOO 



= n''S e-"''^' (st: -l/n)ds-n''6 e""'*'* { - Z/n} rfs 

by the argument above. Since the absolute value of the fractional part does not 
exceed one, we have the estimate 

poo -1 / 1 \ 



j=0 
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This implies 

TTL 

!{z)^-Y^e n r 1 + - z-i^+0(l) (62) 



n ^-^ \ b 

1=1 



m 



Therefore f{z) either blows up like z^t (when X];=i ^ 7^ 0) or it is bounded 

TO 

-27ri — i*" 

(when X;r=i e " =0). 

The Fourier expansion of the fractional part can be used to calculate the transseries 
as we did for /3 = 0, but we shall omit the calculation here. 

For special values of 6, asymptotic information is relatively easy to obtain on a 

r + 1 

dense set along the barrier. This is the case when b — where r S N; then, the 

r 

transseries contains exponential sums in terms of integer powers, fc'', a consequence 

of the duality relation - H — - = 1, which at the transseries level is of the form 
a 

^g-fc 2 _|_ ^ e'^'^ ^ 52 where (71,52 are power series. We illustrate this for 

b^l. 
2 

Without loss of generality, we assume (3 < Q. The transseries of / is given in 
(|87p. To estimate the asymptotic behavior of f{z) as z approaches the imaginary 
line, we rewrite ((87)) as 




3i 105i5 , 
s ^^s- 



00 



4%/2(7ri)5 87r3 256V27ri 
32i7r3fc3 

4V2i~27r , 1 

z 



■ ds 



fe=i 



TT — ' \ Z 







8\/27r^ 



Watson's Lemma implies that 

32i7r3fc3 

"^TT .r-^ 1 

3z 



/W-^^^^^E^^e 27z2 (64) 



fc=i 



The is sum in ([64]) is similar to the sum with 6 — 3, and can be estimated in a 
similar way: 



20 O. COSTIN, M. HUANG 



Setting = V + 2ni— , we have z = -—f=J hi — y + o{y). The 

2"z^ n 3\/3 y m \mJ 

4:7rz / ?T- 

asymptotic behavior can be obtained for (3 = — (this includes all ra- 

3v3 V 



tionals) by substituting y = ztti — m the above estimates. Setting 

27z^ n 



,(5((5>0),a direct calculation shows that 

3^3 Vm 



in "^3 

VM^)m -^(-Y-T. + oil) (65) 



1=1 



3.8. Details of the proof of Theorem [TJ Consider, more generally, 



k=0 



where g{k) > is a real function andj^ \a{t)\'^dt — 
We find the behavior of 



\f{S,/3)\^dp 



/3o 



(/3i-/?o)5]|a(fc)| 



2g-2g(fc)5 



fc=0 



I ^ a{k)a{j) ^-(fi(k) + ci(j))S^(g(:i)~g(k))2nil3o I^ig(j)-g{k))27rt{0i~-l3o) _ ^ 

2^* - ff(fc) 



3o 



(66) 



where /3o,i G K are arbitrary, or after m integrations, 



F{S)^ / ••• / / \fiS,p)fdpdPo---d^„ 



Coci.--C„^|a(fc)pe-^^W^ + Ol^ 



fc=0 



a{k)a{j) 



igU) - 9{k)y 



(67) 



NATURAL BOUNDARIES 



21 



Note that 

a{k)a{j) 



E 



-(g(k)+gU))S 



oc oo 

= 2EE 



fe=0 n=l 



a{k)a{k + n) 



{g{k + n)-g{k)y 

u{k + n) 



-{g{k)+g{k+n))S 



fe=0 



under our assumption. 
If furthermore we have 

oo 

E 

n=l 

then we obtain 



n=l 



{g{k + n)-g{k)) 

OO 



n=l 



_(g(fe+„)-g(fc))5 



a(A; + n) 



{g{k + n)-g{k)y 



(68) 



a(fc + n) 



{g{k + n)-g{k)y 



o{a{k)) 



F{6) = coci ■■■CmJ2 l«(fc)l'e-^«('=)^ + O E 



fc=0 



KM! 



-2ff(fc)5 



where > 1, G{k) ^ oo as fc ^ oo, or 



■2gik)S 



Co • • • Cm + o(l) as ^ — * 0"'" 



(69) 



(70) 



The result now follows from the following lemma. 

Proposition 8. Assume 

(i) hn - [0,oo) are locally L^, and 



(a) lim / hn{x\-\ \-XN)dx\- ■ ■ dxN = 'meas{B) for any box B =\\[ai,bi 

Jb ~.i ■ 

Then hn ^ 1 in the dual ofC[a,p] for any [a, (3]. 



Proof. We first take N = 2, the general case will follow by induction on N. 

Consider the rectangle Be = {a,b — c) x (0,6), < c < b — a. By changing 
coordinates to x + y = s,y = y' , we get that 

/ hndydx = / hn{s)Ta^b;c{^)ds meas{Tafi;c)\ as n ^ oo (71) 

J B J a 

where Tq^(,;c(0 is the function having Tafi-c as a graph, To_(,;c being an isosceles 
trapezoid with lower base the interval (a, b) and upper base of length 6 — a — 2c at 
height 1. 

We also note that the indicator function of [a, 6], lab, satisfies the inequali- 
ties Ta-c,b+c;c > lab > Ta.b-c- Thus, sincc c is arbitrary and /i„ > 0, and both 
meas{Ta-c,b+c;c)i and meas{Ta^b;c) tend to {b — a) as c ^ 0, we have 

rb 

lim / hn{s)ds = {b — a) = meas{[a,b]) (72) 
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In particular, given a < /S, ||^n||Li[a,/3] are uniformly bounded, that is, for some 
C > 1 we have 

snp\\h,,\\Li[^j]<CiP-a) (73) 

n>l 

Since a continuous function on [a, P] is approximated arbitrarily well in sup norm 
by finite linear combinations of indicator functions of intervals, it follows from ()72l) . 
and the triangle inequality that 

K{s)f{s)ds^ V/eC[a,/?] (74) 

J a 

For general N we use (s) = Jg, /i„(s + xi + ... + XN-i)dxi ■ ■ ■ dxN-i and ([7^ to 
reduce the problem to — 1. □ 



The condition 

a{k + n) 



E 



o{a{k)) 



is satisfied, for instance, if 

(1) 3c > so that c < \a{k)\ < c~^, or \a{k)\ decreases to 0. (Note that 
g{k + n) — g(k) = g'{k + in)n, where t 1, and .g'(fc) — > oo as fc ^ oo.) 

(2) 3c, r so that c < a{k) < k^. g'{k) ^ fc^ for some small e > 0. 

4. Proof of Theorem [6] 

In Appendix i jS.SI we list some known facts about iterations of maps. 

Proof of Botcher's theorem, for i28]} . (Note: this line of proof extends to general 
analytic maps.) 

We write V' -^-^ + ^'^zg{z) and obtain 

ff(^) - Igi^') ^\-^ + \^ - 9{-A) + 5(^')] + ^g{z)g{z^) = N{g) (75) 

Let A\ denote the functions analytic in the polydisk Pi g = D x {A : |A| < e}. We 
write (Uni) in the form (seeHH)) 

g = 21N{g) (76) 

This equation is manifestly contractive in the sup norm, in a ball of radius slightly 
larger than 1/2 in Ax, if e is small enough. For A 7^ 0, evidently ij) = is also 
analytic at zero. □ 

Lemma 9. -0 is analytic in Di for all A with |A| < 1. 

Proof. We have 

i;{z) = ^ (x + + ^x/x) F{X); X = ^(z^) (77) 

For small z 0, ip{z) — 0{z) and thus F{ip{z'^)) is well defined and analytic. Note 
that (|77[) provides analytic continuation of tp from Bp2 to Dp, provided nowhere in 
Dp2 do we have ip = —4/ A (certainly the case if p is small). We assume, to get a 
contradiction, that there is a zq, \zo\ = Aq < 1 so that tp{zo) = — 4/A, and we choose 



''Alternatively, and somewhat more compactly, one can prove the result without the intermedi- 
ate steps l|72|l and l)73|l by upper and lower bounding continuous functions by sums of trapezoids. 
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the least Ao with this property. By the previous discussion, ip is analytic in the open 
disk Uyi^. Then we use the "backward" iteration = X^'^'ijj{z)^/{l + ip{z)) to 

calculate ip{z'^) from ^'(•2)1 starting with z = zq. This is in fact equivalent to ()27l) : 
after the substitution x — {—Xy)~^ we return to ([25|) . with a;o = A/4. Using (vi) 
and (vii) of §5.5) it follows that 1/xn 0, that is, ip{zQ ) -f-> 0. This impossible, 
since is analytic and ■(/;(0) = 0. □ 

Proof of Theorem\Sl (i)- We return to ([28]). Taking a € (0,1) m„ = sup{|?A(z)| : 
l^l < a}/'^ we note that 

TOn+i < ^|A|(m„ + y/ml+Amn/\\\) (78) 

The sequence of m„ is bounded by the sequence of M„, defined by replacing "<" 
with "=" in dill). Since i|A|(a; + ^x^+Ax/\\\) < x ii x > A := |A|/(1 - |A|), we 
have limsup„ M„ < A. By the maximum principle, \^{\^ z)\ < A in D x B. Thus, 
by Cauchy's formula in A we have |V'n(z)| < A for all n and z g D. The radius of 
convergence of (pO|) in A is at least one. By §5.5) (vii), the radius of convergence is 
exactly one. 

Indeed, note first that (a) if 1/' is analytic in D then ?/;' in D, otherwise 
tp'{zi) = would imply 'ip'{z^ ) = in contradiction with ip'i^) = ^- This means 
that if there is a zq, V'(^o) = ^4/A, then ^/zq is a singular point of tjj. 

Secondly, any A of the form 1 + ie with small e correspond to c = 1/4 + l/4e^, 
outside the Mandelbrot set. Thus, in the iteration ^7\ . the initial condition yo = 
— 4/A implies t/„ — s- 0. We can now use the implicit function theorem to suitably 
match yn, once it is small enough, to some value of V' near zero. Indeed, the equation 
yn = Azq" + 0(2:0" ) has 2" solutions. This means that for such a zq, using ([77)1 
to iterate backwards and to determine "0(^0) (noting the parallel to ((27|) '). we have 
■0(zo) — —4/A, and by (a) above, tp cannot be analytic in z in D. 

Formula ([5^ follows by straightforward expansion of (I^S)) and identification of 
powers of A. □ 

Proof of Theorem\^ (ii). The stated type of lacunarity of "ipk follows from p2p by 
induction, noting the discrete convolution structure in A:. □ 

Proof of Theorem\dl (Hi). Continuity of f/jfe in D also follows by induction from 
([32| and the properties of T. By dominated convergence (applied to the discrete 
measure |A|"), for A < 1, ip is continuous in D and the Fourier series converges 
pointwise in 9ID). 

To show convergence of the Fourier series of H we only need to show info \ > 0. 
Now, ip clearly cannot vanish for any zq = D, otherwise ipizQ ) = 0, would imply by 
analyticity ijj = 0. If minn^ |V'(p)l — ^ would be small enough, then mine 2 lV'(p)l !i 
O(e^) <^ e, contradicting the maximum principle for z/'ip{z). 

The rest of the proof is straightforward calculation, using the analyticity of ip. □ 

The extension of the small A analysis to higher order polynomials is also straight- 
forward. 

Note 7. The transseries of the Botcher map at binary rational numbers can be 
calculated rather explicitly. This is beyond the scope here, and will be the subject of 
a different paper. A less explicit expression has been obtained in |5j. We note that 
the constant log2{2TT) in (1-17) 0/ [5_ should be (27r)/log2. 
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5. Appendix 

5.1. Proof of Lemma [4l For every to, we write g{s) — gm + o(s^™) where gm is 
a finite sum, an initial sum in the asymptotic series of g. It is straightforward to 
show that g^i{y) has an asymptotic power series as y — *■ oo. Then, in the equation 
g[s) = y we write s = + e where g,n{sm) = V- Then, y = g{s) = g{sm + e) = 

ffm(Sm)+.g'(6e+(.9(Sm)-5m(s„i)) imphcS £ (gm (Sm ) " ^(Sm))/^' (0 = 0(3"™"") 

where g'{£,) ^ ax^. Now 9 is fixed and m is arbitrary, and then the rcsuh follows. 

5.2. Proof of Lemma [71 We have 

/>oo 

y_e-^9(o) =_Vfc(e-^s('^+i)- e-^f = a; / e-^9('^)5'(s)[sjds 
fc=i 

nOO nOO nOO 

= x e~^"[5~i(u)Jd« = a; / e-'"'' g~\u)du - x e-^"{<7"i(u)}dM (79) 

Jg(0) Jg{0) Jg(0) 

and we also have 

pOO POO POO 

X e-"'g-\u)du = - (e~^»("))WM= / e^^^^^^dw (80) 

Jg(0) Jo Jo 



whereas 

/>C30 />oo 

0< / e-^"x5'(u){u}dM < / e-^"a;5'(M)du = e^f(°) 
5.3. Proof of Lemma [Tl By the Fourier coefficients formula we have 



Co= / fiay)dy+ / ydy- / G(a^)dy 
Jo Jo Jo 

Ck= f f{ay)e-'^^'"ydy+ f ye'^'^'^'y dy ~ f G{ay)e~'^^'''y dy 
Jo Jo Jo 

2/c7r log a \ \ log a J \ log a 

I f ("D"^^ -f (-1)"^^ I 1 rf (MO) 

^-^ n!(2fc7ri — nloga) n!(2fc7r« — nloga) log a \ log a' ' 



(81) 



1 2fcn-jp 

Note that since ' {^k^T^^n^i^) = (n ^ p) we have 

(-1)"+! /-""y^ (-l)"+le-"P 1 ^f-;''^^ 

^ n\(2kTTi — nlosa) In ^ n\ loga 

n=0 ^ to ; to 

J_ / —p 2fc7rip 



loga 







e'h-^^-^dt=- r(-- )-r(--— ,1) 



loga Jo loga \ loga loga 
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The above procedure is justified for k in the upper half plane. By analytic 
continuation the expression holds for k real as well. Eq. ((50|) follows. 
We can further resum the series in the above expression by noting that 



2kTri 



log a 



2kT,i -^—^ — log a 
2km 

k=l 



E 



k— — OD 

-1 



2kTTi 



oo oo 

E 

fc=l 



- loga ( z\ 



2kTTi 



(!)■ 



k—~oo 



— log a f z 
2k'Ki 



e '^dt I t '"s" e ^z^'B^dt 
Jo 

\ogJl-[^^^^\e-^dt 



loga 
27ri 



iog«(i-(Tr"' 



2711 

log az 
27ri 
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logK (-sT^ j e-^^ds (82) 



which can be justified by analytic continuation, for the last expression and the sum 
are both analytic, and equal to each other on the real line. The logarithm log^ is 
defined with a branch cut along R~ . 

We finally obtain the integral representation ([5T|) . valid for z in the right half 
plane. 



Remark 8. Since log^ ^— si°s<'^ — iarg 



27r log s 
loga 



2Txi 



riog 

I log 



we actually recover the last term of 
5.4. Direct calculations for & e N integer; the cases b = 3,b ~ 3/2. 
Proposition 10. If b is an integer, the behavior of 

oo 

f{S + 2Trip) = J2 e-^-''(*+2"*'') (Re((5) > 0) 

A;=l 



where (3 = 2'Kim/n, m and n being integers, as 5 approaches is 



f{S + 2mP)^ 



-E' 

1=1 



m 

'2-Ki — 1° 



n r( 1 + - 



5--^ +0{l) 



Therefore f{S) either blows up like S ^ or is bounded. 

r + 1 

The more general case b = where r is an integer can be treated similarly. 

r 

3 4:7rz / n 

In particular, if b — —, for j3 = — ^^^ w — (this includes all rational numbers), we 

have 

oo rp, — : n ,3 i 

^E^ " 
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Proof. In general, to find the asymptotic behavior of f(z) we analyze the functions 

fk{z)= / e-f^-^'^-P^dp, fceZ (83) 



Letting p = q (^-^^ we have 



Next we let s 



h{q) = g + 27rzgiand ^ (f )^ /^^ e-(4)^ /T^X^*^" 

where the contour Ci is a curve from the origin to oo in the first quadrant. 

We can find the asymptotic behavior of fk{z) using Watson's Lemma P|. By 

iterating the contractive map q I — — — I near 0, we can easily see that ^ — 

\ ^ITTj J s 

is analytic in s''^^, which implies — - — tttt is analytic in s**^^ near with no 

h'{h-^{s)) 

constant term. 

Now let's consider the examples b — 3 and b — 3/2; for 6 = 3 we have s ~ h{q) = 

„ . J. . ds 3? „ 15 / 21i fi „, , 
q + Zmqs and — =-:-^s + „ , „ s — 7777— prs • • • . thus, the asymptotic power 
dq Stt-^ 647r° 1287ry 

series is 



■'^ ' \3j 2 120 792 

4 

The branch point of h^^ is located at sq — tt^ — which is between the 

3v 3 

contour Ci defined above and the cc-axis. As we start rotating z from z > 0, we 
have, cf (f8^ . 



j h'{h'Hs)f'^\z) J h'{h'Hs + so)r 

C2 



ds + 2 - e-^°(^)' / TTTj-r, ^ds (84) 



z; Jo h'{h-^is)) \z) Jo /i'(/i-i(s + so)) 

where the contour C2 starts at 00, goes clockwise around the origin, then ends at 
00. Since now 



ds (—Tri) -1, ,1 5 5 f i . , 
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(85) 



Therefore the transscrics is 

f(z) = T(-)z-i - I - — + — + 
2 120 792 



fc=i 



— I k 'tZ * + 



5 J4 , _7 1 
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(86) 



The calculation for 6 = - is similar: in this case ^ = — — , 

2 rf? 4y2(7ri)i 
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^ 87r3 
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256\^7ri 



+ ■ 



and the asymptotic power series is 



f{z)^T(-]z--^-- 



3C( 



1 2 3150^) 
167r2 ^ 240^ 20487r5 



z^ + - 



The exponential sum is slightly different than in the previous case, for now the 
32 

branch point sq = —77=^77^ lies in the lower half plane, which means the contour C2 

can be deformed to [0,+oo) without passing through any singularity. 

We collect the contribution from the branch point only when arg z decreases to 



from 0. Since 

4 

ds 

dq 



-4-^21 2 TT 2 



(s - So) 



I 2 

3 - ■ 3 ■ g^^f 

we have for the exponential part of the sum 
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16\/27r^ 
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Therefore, for small z in the right half plane, the transseries is given by 



fiz) = T - z-^-- 



1 Kil) 
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1 315C(f) 3 
240 20487r5 
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e 27«2 
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4\/2i~5 7r 
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_^^k-iz^ + ...); _| ^ arg^ < -f 



(-A;)52;"^ + 



16%/27r4 



; 5 ^argz^ f 



(87) 



The effect of the exponential part of the transseries affects the leading order when 
^ — > nearly tangentially to the imaginary line. 
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3 

For example, to see the effect of tlie first exponential term for 6 = — , we let 

z 3 ~ e 27ji! z and z = —ire near the negative imaginary line. 

The critical curve along which the power term and the exponential term are of 
9 1 

equal order is 9 ^ „ r^ log -. □ 

g47]-J J- 




Figure 5. The Madelbrot set (drawn with xaos 3.1 [18]). 



5.5. Notes about iterations of maps. For the following, see e.g., [3l Ffl [T6]. 

(i) For |A| < 1, three types of behavior are possible for the solution of (P5|) : if the 
initial condition xq G the connected component of the origin in the Fatou 
set, then a;„ ^ as n ^ oo. Clearly, xq € T \i xq is small enough, li xq ^ J- , 
the connected component of infinity in the Fatou set, then |x„| — > oo. Clearly, 
xq ^ T if it is large enough. Finally, for xq £ OF = J, the Juha set. 
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a (connected) curve of nontrivial Hausdorff dimension invariant under the 
map, {xn}n are dense in J and the evolution is chaotic. 

(ii) J is the closure of the set of repelling periodic points. 

(iii) For polynomial maps, and more generally, for entire maps, J is the boundary 
of the set of points which converge to infinity under iteration. 

(iv) If the maximal disk of analyticity of ip is the unit disk Di, then i/j maps Di 
biholomorphically onto the immediate basin of zero. If on the contrary 
the maximal disk is D^, r < 1, then there is at least one other critical point 
in Ao, lying in V'(9Dr) ~ Jy, the Juha set oflTfl 

(v) If r 1, it follows that ip{dBi) = Jy. 

(vi) By the change of variable Xn = — + 1/2, ([25|) is brought to the "c 
form" tn+i = tf^ + c, c = A/2 — A^/4. The Mandelbrot set is defined as (see 
e.g. [21) 

M ^ {c:tn bounded if to = 0} (88) 

If c e A4, then clearly j/„ in ([^7]) are bounded away from zero. Note that 
<o = corresponds to xq = 1/2 implying xi = —A/4. 

(vii) is a compact set; it coincides with the set of c for which J is connected. 
The cardioid H ^ {(2e** - e2")/4 : t e [0,27r)} is contained in M; see 
[7]. This means {A : |A| < 1} corresponds to the interior of Ai. We have 
\X\ = l^cGdMcM. 

5.6. Overview of Borel summability and transseries. There is a vast litera- 
ture on transseries, Borel summability, and resurgence, see, for example [T9l. Most 
of the modern theory originates in Ecalle's work [3]. 

Definition 11. We say that f is given by a Borel summable transseries for x > v, 
if there exists a /3 g C, a sequence Ck, with Recfc > Ck for some C > 0, and a 
sequence of functions Yk, analytic in a neighborhood o/(0, oo), having convergent 
Puiseux series at zero, and \Yk{p)\ < \B^e'^P\ (where B and v are independent of 
k ) such that 

oo 

f{x) = e-'^'^x^^+^'^I^CYk (89) 

fe=0 

where C is the usual Laplace transform: 

/•oo 

{CY){x) = / e-P^Y{p)dp (90) 

The definition for other directions 9 in the x complex domain is obtained by chang- 
ing the variable to x' = xe^^^ . 

The Borel- Laplace summation operator is denoted by CB. 

Definition 12. A formal power series in powers of 1/x is Borel summable as 
X oo if it is the asymptotic series of CY , where Y is as in Definition llll (We 
note that by Watson's Lemma "21, {CY){x) has an asymptotic series as x ^ oo, 
which is the termwise Laplace transform of the Taylor series of Y at zero.) 

Transseries representations contain therefore manifest asymptotic information. 

Definition 13. A function Y{p) is resurgent in p in the sense of Ecalle 9 , if it 
is analytic on the Riemann surface of C \ J , where J is a discrete set, and has 
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uniform exponential bounds along any direction towards infinity cf. [ 19^ [j By abuse 
of language, f{x) is called resurgent if it satisfies the requirements in Detinition \ll\ 
and all are resurgent. 

This is especially useful when global information about / for x G C is needed: de- 
formation of contours in p, and collecting residues when/if singularities are crossed, 
provides a straightforward way to obtain this information. 
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It is also required that certain relations, called bridge equations, hold; in our case it would be 
easy to derive them from the explicit form of H, but we omit this calculation. 



